Jamming, Yielding and Rheology of Weakly Vibrated Granular Media 
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We establish that the rheological curve of dry granular media is non-monotonic, both in the 
presence and absence of external mechanical agitations. In the presence of weak vibrations, the 
non-monotonic flow curves govern a hysteretic transition between slow but steady and fast, inertial 
flows. In the absence of vibrations, the non-monotonic flow curve governs the yielding behavior 
of granular media. Finally, we show that non-monotonic flow curves can be seen in at least two 
different flow geometries and for several granular materials. 

PACS numbers: 45.70.Cc,81.05.Rm,45.70.-n 



Granular media are collections of macroscopic and 
athermal grains which interact through dissipative, fric- 
tional contact forces and that jam in metastable con- 
figurations Under sufficient shear stress, or when 
mechanically agitated, granular media yield and flow [lol - 
flil ]. Here we ask: what is the rheological scenario that 
connects these observations? 

The classical example of a inclined layer of sand illus- 
trates that granular media exhibit a finite flow thresh- 
old, and once the material yields, the flow rate jumps to 
a finite value [l], 0-11] • The simplest flow scenario that 
captures this behavior is sketched in Fig. la, where in 
analogy to static and dynamic friction, the static yield 
stress a s is assumed to be larger than the dynamic yield 
stress ad- 

Such scenario implies that the flow rate continuously 
decreases to zero when the stress (inclination angle) is 
lowered, as can be seen from following the flow curve 
in Fig. la. In experiments, however, the flow is found 
to stop discontinuously: stress-controlled granular flows 
have a minimal flow rate This calls for a more 

elaborate flow curve than the one in Fig. la. We sketch 
one that was proposed earlier [(| [j| in Fig. lb, where the 
negative slope, signalling an instability, leads to a "for- 
bidden" range of flow rates — when the stress is lowered 
below a m m, the flow rate jumps to zero 0, 

In this Letter, we firmly establish the existence of non- 
monotonic flow curves for granular media by probing 
their rheology, both in the presence and absence of exter- 
nally supplied vibrations of strength T. In the absence 
of vibrations (T = 0), the flow curves indeed are of the 
form as sketched in Fig. lb, i.e., with a finite yield stress, 
and a dip at intermediate flow rates. As we will see, the 
rheology at T — can be seen as a limiting case of the 
more general rheology at finite vibration strength. For 
r > we find flow curves as sketched in Fig. lc, i.e., with- 
out a clear yield stress and with two competing branches 
of positive slope. This lower branch is consistent with 
recent experiments which have found that external agi- 
tations decrease or even quash the yield stress of granular 
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FIG. 1: Hypothetical flow curves for granular media, relating 
strain rate j and stress a. (a) Monotonic flow curve with a 
different static (a s ) and dynamic (ad) yield stress, (b) Non- 
monotonic flow curve - the region with negative slope signals 
'forbidden' strain rates in stress controlled experiments, (c) 
A typical flow curve observed in our experiments for finite 
vibration strengths. 



materials 1(J- 16 1 , allowing for very slow flows under con- 
stant, small, stresses and finite vibration strength. Non- 
monotonic flow curves such as depicted in Fig. lc are 
well-known to arise for polymer melts, micelles and vis- 
cous suspensions [1, IrH - flfjj l but have not been observed, 
to the best of our knowledge, for dry granular media. 

The non-monotonic flow-curves are obtained in experi- 
ments with controlled flow-rates, and negative slopes are 
expected to cause instabilities and hysteresis when the 
stress is controlled [lj| . For < T < 1 , we find that stress 
sweeps through the unstable regime lead to a concomi- 
tant hysteretic transition of the flow rate between the two 
stable branches. For r = 0, we find that stress-controlled 
yielding, i.e., the hysteretic transition from zero flow to 
finite flow rate, is intimately connected to the dip in the 
r = flow-curve. We finally show that non-monotonic 
flow curves like in Fig. lb can be seen independent of 
flow geometry and for several, but not all, granular ma- 
terials. Our findings of the robust connection between 
non-monotonic flow curves and rheological instabilities 
shine new light on the nature of the granular jamming 
and yielding transitions. 

Setup — Our experiments are carried out in a split- 
bottom shear cell (Fig. [2]), in which a layer of glass beads 
(diameter 1-1.3 mm) of depth H = 18 mm is driven 
by the rotation of a rough disk of radius R s — 4 cm 



2 




10" 

Q [rps] 



container 



FIG. 2: (Color online) Flow curves T(fi) for T = 0, 0.1, ... 1 
as indicated. We normalize T by the dynamic yield torque 
Td- Inset: vibrated split-bottom rheological setup. 



mounted flush with the bottom. This flow geometry has 
been studied extensively and produces smooth, robust 
and well-controlled granular flows (20l - [22l |. 

An important novel aspect of our setup is that the 
shear cell can be vibrated vertically and sinusoidally (dis- 
tortion < 1%) with frequency / and amplitude A. We fix 
/ = 63 Hz in the middle of a frequency window where no 
mechanical resonances arise. Strictly vertical vibration is 
ensured by guiding the motion of the shear cell with a lev- 
elled square air-bearing (4"x4", New Way) which is cou- 
pled to a electromagnetic shaker (VTS systems VG100). 
We control the dimensionless shaking strength, defined as 
T := A(2irf) 2 /g, with a feedback loop to within < 10~ 3 . 

We also control the rotation rate, f2, and applied 
torque, T, by a rheometer (Anton Paar DSR 301), which 
is coupled to the vibrating cell by means of a flexure 
with a torsional spring constant of 4 Nm/rad and com- 
pressional spring constant of 5 x 10 2 N/m. We perform 
rheological experiments at fixed T and either control the 
torque T and measure the resulting rotation rate or 
vice-versa. All flow experiments are preceded by ap- 
propriate pre-shear. Disk rotation is always continuous; 
stick-slip is not observed. Note that, as in other flow 
geometries, the local strain rate and stress in the split 
bottom cell vary throughout the cell [23| . We thus probe 
the grain rheology with T as a proxy for the stress a, and 
f2 as a proxy for the strain rate [7|, |22|. Hence, the ex- 
perimentally observed curves for T(fi) are best thought 
of as global flow curves. 

Main Phenomenology: Flow curves — Fig. 2a shows 
the flow curves T(f2), determined in experiments in which 
the rotation rate is controlled, and the average torque 
T is measured (after removing transients). 

The flow curve for T = is non-monotonic. For small 
flow rates (fi < 10~ 3 rad/s), the stress reaches a plateau 
from which we determine the dynamical yield torque Td 
as 13. 9 ±0.1 mNm — this value is set by the geometry [2(| 



and the effective friction coefficient of the grains. For 
increasing fi, T decreases until it reaches a minimum 
torque T m i n of about 12.1 ±0.1 mNm at f2 ~ 0.3 rps. This 
non-monotonic effect is substantial in magnitude and has 
not been observed for granular flows before. Around the 
minimum, the inertial number near the split is of order 
one, and we associate the increase of torque for larger 
rates with the onset of inertial flows 0, H| ■ 

The flow curves for T > exhibit similarly non- 
monotonic behavior, but differ for small f2. As shown 
in Fig. 2a, the flow rates over which the flow curves have 
negative slope become smaller for larger T. At the lower 
Q range of this regime, T(f2) reaches a local maximum 
and for even smaller f2 < 0.02 rps, we observe a decrease 
of T with fl as T ~ log(fl). In additional experiments at 
fixed T we have carefully checked that the flow is stable 
and steady in this positive slope regime. This regime only 
exists for finite agitation strength and signals a novel flow 
regime of mechanically agitated granular flows which is 
unique to T > 0, as suggested in Fig. lc. 

We conclude that our flow curves are consistent with 
the flow scenarios depicted in Fig. lb-c. In other systems 
with non-monotonic flow curves, fixing the flow rate in 
the negative-slope regime typically leads to a separation 
of the system into two regimes, one with low, and one 
with large strain-rate: in other words, shear banding fl9| . 
In contrast, we have not seen any clear evidence for such 
behavior in our system — the flow profiles as observed 
at the free surface do not appear to change when we fix 
the flow rate in the negative-slope regime. We note that 
the standard shear banding mechanism depends on the 
shear stresses being sufficiently homogeneous, while in 
our system we have a strongly inhomogeneous stress field 
emanating from the split in the bottom (23j — this inho- 
mogeneity is crucial in obtaining a smooth granular flow, 
but may hinder the observation of additional shear band- 
ing. Our set up does not allow us to determine whether 
dilatancy plays a significant role in the development of 
the non-monotonic flow curves. Yet, since we only con- 
sider slow steady flows, the usual inertial and transient 
dilatancy effects 0, HiJ are certainly ruled out. 

Rheological instability for T > — We now turn 
our attention to torque-controlled experiments, and will 
probe whether the negative sloped regime of the flow 
curves for T > leads to hysteresis. To do so, we 
slowly ramp T /Td up and down between 0.8 and 1.1, i.e., 
through the multi-valued regime. Fig. |3k illustrates the 
resulting hysteresis loops. Ramping upwards, we observe 
a sudden jump from the slow, mechanically agitated flow 
branch to the rapid, inertial branch. Ramping down- 
wards makes the flow rate jump back to the slow, me- 
chanically agitated branch — there is considerable hys- 
teresis between the stresses where these jumps happens. 
For smaller T, the gap between slow and rapid flow rates 
increases, consistent with the flow curves shown in Fig. 2. 

In Fig. [3Jd we further strengthen the direct connection 
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FIG. 3: (Color online) (a) Finite V hysteresis loops for T — 
0.1,0.3,0.5. (b) Several hysteresis loops at T = 0.5, with the 
rheological data from T — 0.5 from Fig. 2 overplotted in black. 



between the negative slope of the T(f2) curve and the 
hysteresis observed in the ft(T) curves for the example 
of T = 0.5. We combine several torque controlled data 
sets with the appropriate flow curve, and observe that 
while the precise location of individual hysteresis loops 
fluctuates, the characteristic torques remain confined to 
an interval which coincides well with the minimum and 
maximum of the T(f2) curves. We concluded that for 
r > 0, hysteresis and negatively sloped flow curves are 
directly related. 

Static yield vs dynamic yield for T — — In the ab- 
sence of vibrations, the connection between flow curves 
and instabilities is more subtle, as there is only one sta- 
ble branch with a finite flow rate (although one could see 
the jammed state at small stress as the second "stable 
branch"). We now ask the following question: when we 
slowly ramp up the torque, how is the resulting yielding 
behavior influenced by the non-monotonic flow curve? 

To answer, we study the statistics of static yielding by 
ramping up the applied torque at 0.5 mNm/s, fixing T = 
0, and measuring the ensuing il(t). We identify a yield 
event whenever > Q,Th — 0.016 rps. Our statistics are 
robust for VIth between 0.002 and 0.13 rps, and over the 
duration of the experiment (~ 10h), we do not observe 
any appreciable drift in the properties of the distribution. 
As indicated in Fig. H^,, we observe two types of yielding: 
micro yielding, where VL only briefly peaks above flxh 
(diamond), and global failure, where the increase in is 
dramatic and persistent (square). 

We measure the statistical properties of both types of 
yielding over 1943 torque ramps. In Fig. 0t> we show the 
probability distribution functions of the micro yielding 
torques (gray) and the global yielding torques (black). 
Consistent with the flow curve for T = (Fig. 3c), we ob- 
serve that micro yield events do not occur above T/Td = 1 
— once the torque is above the dynamic yield threshold, 
the material flows at a rate given by the positive slope 
region of the flow curve. What is surprising is that global 
yielding can happen both below and above Td — invali- 
dating the simple picture of a lower dynamic, and higher 
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FIG. 4: (a) Typical examples of 0.(T) curves during a torque 
ramp. The dashed line indicates the threshold rotation rate 
^Th = 0.016 rps that defines a yield event, o indicates a 
micro yielding event, □ a global failure event, (b) The proba- 
bility distribution function of micro yielding (grey) and global 
yielding (black). Grey bars indicate the onset value for the 
rheological instability T m i n , and the dynamic yield stress Ta- 
lc) Flow curve for T = 0. 



static yield threshold. Rather, the probability for com- 
plete failure is bounded by T = T m i n — for lower torques 
the flow curve shows there is no steady state flow, and 
only micro yielding events are observed. 

We conclude that, due to the non-monotonic flow 
curve, there are three yielding regimes. Below T min , mi- 
cro yielding without global yielding; above T min but be- 
low Td, micro yielding with a finite probability of global 
yielding; above T = Td all yield events become global. 
This intermediate regime is clearly inconsistent with the 
simple picture where hysteresis in jamming and yielding 
of granular media is explained from the static threshold 
exceeding the dynamic one (Fig. la). 

Robustness non-monotonic flow curves — How specific 
to the precise flow geometry and materials used are the 
non-monotonic flow curves we observe? We test this in 
two ways. First, we can test the dependence of the rhe- 
ology on the shape of the shear bands in the split bot- 
tom cell itself with the glass beads mentioned earlier: for 
larger filling height, i.e. H/R s = 0.7, the flow field is 
qualitatively very different 23-13]- We have conducted 
a series of experiments and we observe qualitatively the 
same flow curves and instabilities (not shown here). 

Second, we probed the existence of the non-monotonic 
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glass beads 3 mm 
plastic flakes 1-1.5 mm 
aqua collisa 0.5-2 mm 
bronze flakes 0.5 mm 
pmma beads 3 mm 
aluminum oxide- 
flakes 0.5-3 n 
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FIG. 5: Rheological curves for a range of materials as indi- 
cated, where a vane (see sketch) directly coupled to an Anton 
Paar MCR501 rheometer was placed just above the bottom. 
All torques are normalized by their mean value and (with the 
exception of the bottom curve) shifted upwards for clarity. 



flow curves for a range of materials in a standard vane 
geometry at T — 0. The top curve shown in Fig. [5] shows 
that the non-monotonic flow curves for glass beads at 
r = are not specific to the split-bottom geometry. How- 
ever, details of the flow curves are material dependent - 
some materials (aluminum flakes and PMMA beads) do 
not have a negative slope at all, while others (bronze 
flakes, aqua collisa) have the positive/negative/positive 
slope combination that we saw above for V > only. We 
do not know the cause of this material dependence, but 
suggest that material dependent plastic flow in the con- 
tact asperities may destroy the negative slope, and that 
for such materials, temperature may play a similar role 
as vibrations does for glass beads. 

Conclusion — We have shown that non-monotonic flow 
curves are a robust feature of slow granular flows. As a 
consequence, the yielding transition at T = exhibits 
all the hallmarks of a 1st order or subcritical transition. 
For finite T, the yield stress vanishes but the hysteretic 
nature of the transition persists up to a large value of 
r (around 1 for this filling height, somewhat smaller for 
larger H/R s ), where the transition becomes 2nd order 
and smooth. 

Our data suggest that material properties are key, and 
while we cannot exclude collective effects 
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a sim- 
ple decrease of the particle-particle dynamic friction co- 
efficient with rate may be sufficient 30]. Our findings 
provide crucial input for refiningand extending current 
descriptions of granular flows 
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